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CN I We analyze the ferromagnetic Ising model on a scale-free tree; the growing random network 
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Mh' model with the linear attachment kernel A^ = k -\- a introduced by [Krapivsky et al.: Phys. 
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Rev. Lett. 85 (2000) 4629-4632]. We derive an estimate of the divergent temperature T^ below 
which the zero-field susceptibility of the system diverges. Our result shows that Tg is related to 



1-^ ' a as tanh(J/Ts) = a/[2{a + 1)], where J is the ferromagnetic interaction. An analysis of exactly 



solvable limit for the model and numerical calculation support the validity of this estimate. 
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I. INTRODUCTION 



Many real-world systems range from the structure of Internet or WWW to social rela- 
tionship between human society, or prey-predator relationship in food webs are described 
topologically as scale- free networks (SFNs) [l|, l2, Isj. In a SFN, the degree distribution 
P{k), where degree k is the number of edges connected to a node, has a power-law decay 
P{k) oc k~^. The degree exponent 7 takes 2 < 7 < 3 in many real networks |^. SFN studies 
have been carried out actively in recent years, including various processes taking thereon, 
e.g., network failure, spread of infections, or interacting systems, which attract numerous 
applications and fundamental interests about critical phenomena [2|, y, |5|. Among them, 
the ferromagnetic Ising model on SFNs is a basic prototype to understand how network 
topology influences cooperative behaviors on SFNs. Indeed, previous studies have demon 
strated the extreme influence of the network topology by both analytical [g, |7|, |8|, |9|, 
and numerical methods 



. _ „ „ . lo|,lii| 

[12I, |l3|, [ij] . For example, Dorogovtsev et al. [9\ analyzed the Ising 



model on an uncorrelated SFN model with a power-law degree distribution by the Bethe 
approach to show that the phase transition exists at a finite temperature for 7 > 3, while 
the system remains in the ferromagnetic phase at any finite temperature for 7 < 3, and its 
critical exponents vary depending on 7. 

In the previous approaches, rather dense part of a SFN has been analyzed in details, 
while our knowledge of how the spins on leaves (nodes whose degree is one) behave is still 
missing. In practice, many real SFNs have a number of leaves. How does such network 
topology influence the critical behaviors? This paper focuses on a simple case; the Ising 
model on a tree-like SFN which includes many leaves, but has no cyclic path. We note that 
the ferromagnetic Ising model on trees behaves quite differently from that on regular lattices; 
the spin system on a tree has no spontaneous magnetization at any finite temperature, while 
its zero-field susceptibility remains to diver ge b elow a certain temperature Tg (we call it the 



divergent temperature) 
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22l |. In other words, the leaf-spins are 



extremely sensitive to the external field below T,. Our aim is to clarify the relation between 
the divergent temperature T^ and the degree distribution P{k). 

In this paper, we analyze the Ising model on the growing random network (GN) model 
introduced by Krapivsky and coworkers 23j. The GN model probabilistically generates a 
sampled tree Tn with A^ nodes as follows. One starts with one isolated node. At each time 



step, a new node is added and links to a preexisting node. The probability that the new node 
links to a node with degree k is proportional to the attachment kernel A^. In this paper, we 
focus on the linear attachment kernel A}, = k + a {a > —1). The degree distribution of a 
resulting tree is determined by the offset a in the attachment kernel. For the infinite offset 
a -^ oo, the degree distribution reduces to the exponential form P{k) = 2^^ . For a finite 
offset a, the degree distribution satisfies a power-law decay P{k) oc A;"''', where the degree 



exponent is 7 = 3 + a [23|. As the offset a is smaller, the resulting degree distribution 



P{k) is more fat-tailed. Particularly, the GN model with the no offset a = corresponds 



to so-called Barabasi- Albert tree 
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27|, which has the degree exponent 7 = 8. We 



investigate how the divergent temperature Tg is related to a to show the extreme sensitivity 
of trees with fat-tailed degree distribution (2 < 7 < 3). 

This paper is organized as follows. In section [TTl we introduce the Ising model on the GN 
model. In section llllt we derive the exact expression for the system susceptibility and the 
divergent temperature for the infinite offset case a -^ 00. In section HVl we give an estimate 
of Ts for the general offset case. In section |Vl we show our numerical results to support the 
validity of our estimate, and suggest that for the no offset case, Tg diverges and an unusual 
scaling exists. Section IVII is devoted to summary. 

II. MODEL 

In this section, we introduce our model; the Ising model on the GN model. The Hamil- 
tonian is as follows: 

n=-JY,S^SJ-hY,S^, (1) 

where J(> 0) is the ferromagnetic interaction, h is the external magnetic field, and Si{= ±1) 
is the Ising spin variable on the node i. The first sum is over all edges of a network, 
and the second one is over all nodes. In the following sections, we calculate the zero-field 
susceptibility of this model. The zero-field susceptibility is expressed in terms of the spin- 
spin correlation functions as Xsys = -^ X^i j=i('S'j5'j), where the angular bracket denotes the 
usual thermal average. For trees, the correlation function between two Ising spins Si and Sj 



on a sampled tree Tat is given as 



19| 



{SiSj)=t''''^^^''\ (2) 



where t = tanh(/3J), P = 1/T, T being the temperature, and dij(TN) is the path length 
between the node i and j on the tree Tn- Accordingly, the one-spin susceptibility XiO^n) 
of a spin on the node i of a sampled tree T^r is 

N 

X.(T^) = /3 5^t'^-(^-\ (3) 

and the system susceptibility 

N ^ N 



Xsys 



™ = ^E^«(T^) = ^E^'"^^^^- (4) 






Note that the system susceptibility is related with the so called average correlation volume 

^V m- ^V =TXsys- 

III. INFINITE OFFSET CASE 

In this section, we consider the GN model with the infinite offset a ^ oo. We derive 
the exact form for the mean system susceptibility x^(A^, T) and the divergent tempera- 
ture Tg. Here A = ^rp P{Tn)A{Tn), and P(Tn) is the normalized probability of a tree 
Tiv- Suppose that Tat+i is created by attaching the (A^ + l)-th node to the n-th node 
of a preexisting tree Tjv. Then the distance from the new node to all others is given as 
di,N+i(T^N+i) = c^Af+i.ilTAf+i) = 1 + (i„^j(Tjv) for 1 < z < A^ [27]. The diagonal element is 
zero: c?jv+i jv+i(TAr+i) = 0. Note that for trees, a path between each two nodes is unique, 
so the matrix elements do not change their values once formed during the growth process. 
So we obtain the recursion relation for the averaged total susceptibility PVn = Nxsjs{N, T) 
as follows: 

N+l 

N N 

= Y ^(Tiv+i) Yl ^'^'■'■^^'^+'^ + 2 XI ^(T^v+i) 5^^'^^'^+'^^'^+'^ + 1- (5) 

Tjv+i ij=l Tjv+i J=l 

The first term of the r.h.s. is 

N N 

Y ^(T^+i) Y ^*'^^^"^^^ = Y ^(T^^) Yl ^''"^^"^ = ^^' (6) 



and the second term of the r.h.s. is 

N N N 



l-iV+1 



i=l Tw n=l 1=1 



Here we use P(Tn+i) = P(TAr)P(n|TAr), where P(n|TAr) is the conditional probabihty that 
the newly-added node links to a preexisting node labeled as n on the tree Ttv- Combining 
Eqs.dnD and ([7]) with Eq.([n]) gives the evolution of Vat as 

Af N 

Vn+1 = VN + l + 2tY, ^(Tjv) Yl P(^\^n) J2 t"'-^^''^- (8) 

For the general offset case, it is hard to solve Vat expicitly since the probability P(r2 |Tjv) 
is proportional to the kernel A^. Fortunately, the infinite offset case is within reaching 
distance. In this case, the conditional probability is independent of which node is attached: 
P(n I Tat) = -^ for any n. Thus Eq.® is evaluated as 

9f ^ ^ / 9f^ 

vw+i = v^^ + 1 + 1 E p(T^) E E ^'^'"^^"^ = 1 + (i + 1) ^^- (9) 

Tiv n=l i=l 

This recursion equation is solved explicitly to obtain the mean system susceptibility as 

^ m=2 ^ ' k=l 

By evaluating the temperature below which the system susceptibility fITOl) diverges, we find 
that the divergent temperature is given as (see appendix [Al) 

tanh( J/T,) = ^. (11) 

Moreover, expanding Eq. llTOl) around the divergent temperature, we obtain the finite size 
scaling form for the infinite offset case: 

Xi^(iV,T) ^ (logiV)/[(T-T,)logiV], (12) 

where f{x) is a scaling function in this case. 

IV. GENERAL OFFSET CASE 

In this section, we proceed to the general offset case. We give an estimate of Tg by 
calculating a lower bound of the system susceptibility. First, we can calculate approximately 



the one-spin susceptibility of the initial node xi- Our calculation is based on a mean field 



approach by Szabo et al. 



261]; the original stochastic model is approximated by a uniform 



branching tree where every node on any level has the same number of offsprings. Let ny 
denote the number of nodes at the l-th level, which means the distance from the initial node 
is /, on the tree with N nodes. When the new node is added to the tree with N nodes, the 
probability that the new node links to any node at the /-th level is 

^N + ^jv + any 
{2 + a)N-l 

Here the new node is stationed at the (/ + l)-th level, so we obtain 

<V.' = n'-»' + ''l-jy^ (/>1), (13) 

where ci = 1 + a, C2 = 2 + a, and the initial condition is n}^ = 1 for all A^. Now we 
introduce the generating function 

oo 

G^(t) = J]4V. (14) 

1=0 

Note that G^i^t) is related to Xi as Xi = PG^it). From Eq. (TT3|) . we obtain the recursion 
relation for the generating function as follows: 

(c2iV - l)GN+i{t) = {C2N + cit)Gjv(t) - 1. (15) 

It is easily solved that 

^"-' + 'r(i + c2-c,t)^^^r(M + i-c2-)- ('^) 

For A^ ^ 1, the summation of the second term can be approximated as 

Er(M + C2"^Cit) v^ ,^£li±i_i ,,£ii±i ,^„s 

r(M+l-C2-i) ^ ' ^ ^ 

M=l ^ ^ ^ Af=l 

SO that 

GN{t)^N^. (18) 

Thus we obtain the one-spin susceptibility of the initial node as 

Ti{N,T) = pG^it) ~ (3N^^^f^, (19) 



which diverges for any T and any a{> —1). 

Now we evaluate a lower bound of the system susceptibility. We consider a subtree which 
consists of a node at the s-th level and its descendents. The number of node at the (s + /)-th 
level is given as n^ I^n ■ Among the total susceptibility of the subtree, the contribution 
from the node-pairs whose paths go through the level s is G^^, — [rv^ /'"'iv ]^^^Af s+i- Here 

M+i) 

N 



Gn, = Y.%v''^ (20) 



1=0 f^N 



corresponds to the one-spin susceptibility of the node at the s-th level. The total suscepti- 
bility of the whole tree is evaluated as 






NTXsys{N,T) — y ^ [G^,. - t Cjv g+i — (T^Ji'-N 

s=0 ^Af 

oo 

= nPG%^, + {l-t')J2ni^'^G%,^^. (21) 

s=0 

The second term is non- negative for < t < 1, so we obtain a lower bound of the system 
susceptibility as 

Tx^s{N, T) > l^n^^G^, = ^G], ^ N'"^-'. (22) 

Note that the exponent includes t. By evaluating where this bound diverges, we obtain an 

estimate of T^ : 

tanh( J/T.) = .rr^y (23) 

2(a + 1) 



which reduces the exact solution ( TTT1) for the infinite offset « — > cxd. This relation indi- 
cates that as the offset is smaller, the divergent temperature is higher. Particularly, we 
immediately find that T, diverges for — 1 < a < 0. 

V. NUMERICAL CALCULATIONS 

In this section, we calculate the zero-field susceptibilities numerically for the GN models 
with several values of offset a. We generate trees for a given offset to calculate the suscepti- 
bilities by using Eq.([2]). First, we show the results for the infinite offset case a — * oo. Figure 
[U compares the numerical result for the system susceptibility with the analytical forms (TTOl) . 
For convenience, we set J = 1. We find that the analytical forms agree well with the nu- 
merical ones. Figure [2] plots the finite size scaling around the divergent temperature. The 



number of nodes is taken from 2^° to 2^^. The system susceptibilities are averaged over 100 
samples. We find that the scaling works quite well. 

Next, we turn to the finite offset case. Figures El plots the finite size scaling around our 
estimate Tg for the mean system susceptibility with the offset 6, 4, and 1. As a result, we 
find that our finite size scaling similar to that for the infinite offset (1121) is quite well fitted 
as long as an offset is not small. These results support that our estimate gives the exact 
divergent temperature. On the other hand, our scaling does not work well for the small 
offsets, e.g., a = 1 or 2, where scaling exponents there may depend on the offset strongly. 

Finally, we consider the no offset case. In FigJl]-(a), we plot the mean system sus- 
ceptibilities log[x^(A^, T)] with several nodes from A^ = 2^^ to 2^^. Now we rescale these 
susceptibilities as \og[Tx^{N, T)]/ log A^. The rescaled system susceptibilities are quite well 
fitted for very wide temperature range as seen in Figjl]-(b). This indicates that Tx^{N,T) 
goes to the infinity as A^ ^ oo. In addition, our result means that the following unusual 
scaling for the system susceptibility (in other words, for the average correlation volume ^y) 
is satisfied: 

Tx^(Ar,T) = [^(T)]'°s^, (24) 

where g{x) is a scaling function. Unfortunately, we have not obtained the derivation of this 
scaling yet. But this relation is derived partially by the following approximation. Bollobas 



and Riordan 2J] derived that the expected number Ei of shortest paths of length / for the 



Barabasi- Albert tree with A^ nodes is given as 



A^^ 1 (i-i°giv)^ , , 

El ~ f 2iogiv (25) 

2 V2vrlogAr ^ ^ 



for A^ ^ 1. Using this distribution, we approximate the system susceptibility as 

N-l 

e 2iogjv fdi 



2 ^/2n\ogN Jo 

for logt > — 1. This approximation shows the system susceptibility Txsys{N, T) holds the 
scaling relation ( 12^ at least in a low temperature region. Interestingly, this scaling form 
remains to be satisfied even at a high temperature where this approximation (126|) is not 
valid. 

8 



VI. SUMMARY 

In this paper, we investigated the zero-field susceptibihty of the Ising model on the GN 
model with the attachment kernel A^. = k + a. Our main finding of this paper is that 
the divergent temperature Tg of the GN model with the offset a is given by tanh( J/T^) = 
a/2{a + 1). The exact expression of the susceptibility for the infinite offset, and the finite 
size scaling of the susceptibilities for the finite offsets support our estimate is exact. The 
finite size scaling form (TT2!) implies that log N can be regarded as the characteristic system 
length L, so that flT2|) can be read as Xs^ ~ LP' f{/S.TL~) with •y = 1, u = 1. 

Our result means that as the offset a is smaller, the divergent temperature Tg is higher 
(Fig|5]). Particularly, T^ diverges when a < 0, that is, the degree exponent 7 < 3. As is 
pointed out in [5], a long-ranged spin correlation covers the whole system below Tg if we 
apply a small local external field on the node i. So one finds that applying a small local 
field, or maybe adding a few shortcuts to spin systems on a tree with a fat-tailed degree 
distribution induces a long-ranged ordering at any finite temperature. Interestingly, our 
result shows that the susceptibilities for the no offset obey unusual scaling IHM . We will 
investigate the origin of this feature in the future. 
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APPENDIX A: THE DERIVATION OF THE DIVERGENCE TEMPERATURE 
OF THE GN MODEL WITH THE INFINITE OFFSET 

In this appendix, we derive the divergence temperature of of the GN model with the 
infinite offset a — > 00. We rewrite the system susceptibility fITOl) as follows: 

N-l ^ 

VN = l+t + 2ty^— —rQm, (Al) 

m=2 ^ ' 

where 

fc=i 



Now, we can derive that the divergence occurs at t = 1/2. At t = 1/2, one find 

Qm=l[ ^ = m, (A3) 

fc=i 

to reduce the system susceptibihty ( lAll) to 

7V-1 ^ N ^ 

VN = l + t + 2ty" - = \"—, (A4) 

m=2 m=l 

i.e., to the harmonic series. Thus, the system susceptibihty diverges at t = 1/2 in the 
hmit A^ ^ oo. We find immediately that the system susceptibihty diverges at least for 
t > 1/2, since vn increases monotonously with t for any A^. Now we show that the system 
susceptibility cannot diverge for t < 1/2. From Eq. (]A2l) . we obtain the following relations 
for Qm'. 

Af-l N-1 

m=2 ^ ' m=2 

and 

2t 

Qm+l ^ Qrn = — Qm- (^^j 

m 
The iterative substitutions of Eqs. flASI) and flA6l) allow one to the following transformation: 

N-l ^ 

^-^ 171(171 + 1 
m=2 ^ ^ 

N-2 



^ ' ^ 2t 1 

-<52 + 2. — 7 TTT'^m - -t:Qn-i 

2 -^^ m(m + 1) A 



m=2 

l + 2t^ \-' (2^)' ^ 2t ^ 1 

2 -^ m m + 1 A/ — 1 A* 

m=2 ^ ' 

....45:p,)-_5:m_q„. (AT) 

m=2 m=2 

In the end, one find 

m=l m,=2 m=l 

This upper bound converges for t < 1/2, so vn doesn't diverge there. Thus, we find that 
the divergent temperature Tg is decided by tanh( J/T^) = 1/2. 
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FIG. 1: The system susceptibility for the infinite offset by the analytical results (jIOp (lines) and 
the numerical ones (points). The number of nodes are taken N = 2^''(red), 2^^ (green), 2^^ (blue), 
and 2^'^ (black), respectively. The average is taken over 100 samples. 
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FIG. 2: Finite size scaling ()12p of the system succeptibility around Ts for the infinite offset. The 
number of nodes are taken N = 2"'^'^ (red), 2"*^^ (green), 2"^^ (blue), and 2"*^^ (black), respectively. The 
average is taken over 100 samples. 
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FIG. 3: Finite size scaling of the system succeptibility around our estimate Tg (|23|) for the offset 
a = 6(top), 4(center), and l(bottom). The number of nodes are taken from N = 2^^ to N = 2^^. 
The average is taken over 100 samples. 
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FIG. 4: (a)The system susceptibility logxi^ and (b) log{Tx^) / log N on the GN model for the 
no offset a = 0. The number of nodes are taken from A^ = 2^*^ to 2^^. The results are averaged 
over 100 samples. 
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FIG. 5: The relation between the divergent temperature Tg and the offset a of the Ising model on 
the GN model. The red-solid line denotes the relation (123p between a and Tg. The blue-dot line 
denotes Tg for the infinite offset a. In the blue-colored region, the system susceptibility diverges. 
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